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Abstract 

In this note we present numerical simulations of binary mixtures. We study the 
diffusion of particles and the response to an external driving force. We find evidence for 
the validity of the Cugliandolo Kurchan off-equilibrium fluctuation dissipation relation. 
The results are in agreement with the predictions of one step replica symmetry breaking 
and the dependance of the breakpoint parameter on the temperature coincides with that 
found in simple generalized spin glass models. 

The behaviour of an Hamiltonian system (with dissipative dynamics) approaching equilib- 
rium is well understood in a mean field approach for infinite range disordered systems ^]. [| . 
In this case we must distinguish an high and a low temperature region. In the low temperature 
phase the correlation and response functions satisfy some simple relations derived by Cuglian- 
dolo and Kurchan [I]. In a previous note j|, |5| we have found some preliminary evidence 
for the validity of these relations. Here we perform a different and more accurate numerical 
experiment and we are able to study the temperature dependance of the phenomenon. 

Let us define our notations. We concentrate our attention on a quantity A(t), which 
depends on the dynamical variables x(t). Later on we will make a precise choice of the 
function A. We suppose that the system starts at time t = from a given initial condition and 
subsequently it is at a given temperature T. If the initial configuration is at equilibrium at a 
temperature T' > T, we observe an off-equilibrium behaviour when we change the temperature. 
In this note we will consider only the case T' >> T (in particular we will study the case 

r = oo). 

We can define a correlation function 



and the response function 



C(t,t w )=<A(t w )A(t + t w )> (1 
SA(t + t w ) 



G(t, ty 



(2) 

e=0 



5e(t w ) 

where we are considering the evolution in presence of a time dependent Hamiltonian in which 
we have added the term 

' dte{t)A{t). (3) 
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The off-equilibrium ffuctuation dissipation relation (OFDR) states some among the corre- 
lation functions and response function in the limit t w going to infinity. The usual equilibrium 
fluctuation dissipation theorem (FDT) tell us that 

G(t) = V^, (4) 

where 

G(t)= lim G(t,t w ), C(t)= lim C(t,t w ). (5) 
It is convenient to define the integrated response: 

R{t,t w )= f drG{r,t w ), R{t)= lim R(t,t w ), (6) 

JO t w ^QO 

which is the response of the system to a field acting for a time t. 
The usual FDT relation is 

R(t)=f3(C(t)-C(0)). (7) 

The off-equilibrium fluctuation dissipation relation state that the response function and the 
correlation functions satisfy the following relations for large t w : 



rC(0,t w ) 

R(t,t w ) « p / X(C)dC. 

JC(t,t w ) 



In other words for large t w if we plot R versus C the data collapse on the same universal curve 
and the slope of that curve is X{0). The function X(C) is system dependent and its form 
tells us interesting information. 

We must distinguish two regions: 

• A short time region where X(C) = 1 (the so called FDT region) and C belongs to the 
interval / (i.e. C\ < C < C 2 .). 

• A large time region (usually t = 0{t w )) where C I and X{C) < 1 (the aging region) 
§0- 

In the simplest non trivial case, i.e. one step replica symmetry breaking [§, |§ , the function 
X{C) is piecewise constant, i.e. 

X(C) = m for C G /, 
X{C) = 1 for C i I. (9) 

In all known cases in which one step replica symmetry holds, the quantity m vanishes 
linearly with the temperature at small temperature. It often happens that m = 1 at T = T c 
and m(t) is roughly linear in the whole temperature range. The relation § has been numerically 
verified in ref. The previous considerations are quite general and can be applied also to 
systems without quenched disorder p"T |. 

If replica symmetry is broken at one step, the value of m does not depend on the observable 
A and the same value of m should be obtained for all the observables. In this case the OFDR 
has an highly predictive power because the value of m may be measured by using quite different 
quantities. 
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The aim of this note is to show that for a binary mixture of spheres the function X(C) 
is similar to predictions the one step formula @ with a linear dependance of m on the 
temperature. It was already shown that in this model simple aging is well satisfied and 
some indications for the validity of the OFDR were already found 0, [| by looking to the 
correlations of the stress energy tensor. In this note we study a different observable and we 
find much more accurate results which confirm the previous findings. 

We consider a mixture of soft particles of different sizes. Half of the particles are of type 
A, half of type B and the interaction among the particle is given by the Hamiltonian: 

where the radius (a) depends on the type of particles. This model has been carefully studied in 
the past |L2], [L3[ [14], |15| . It is known that a choice of the radius such that ub/ua = 1-2 strongly 
inhibits crystallisation and the systems goes into a glassy phase when it is cooled. Using the 
same conventions of the previous investigators we consider particles of average diameter 1, 
more precisely we set 



a\ + 2(a A + a B f + o% 



1. fll^ 
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Due to the simple scaling behaviour of the potential, the thermodynamic quantities depend 
only on the quantity T 4 / p, T and p being respectively the temperature and the density (which 
we take equal to 1). The model has been widely studied especially for this choice of the 
parameters. It is usual to introduce the quantity T = [3 4 . The glass transition is known to 



happen around T c = 1.45 [13 



Our simulation are done using a Monte Carlo algorithm. We start by placing the particles 
at random and we quench the system by putting it at final temperature (i.e. infinite cooling 
rate). Each particle is shifted by a random amount at each step, and the size of the shift is 
fixed by the condition that the average acceptance rate of the proposal change is about .4. 
Particles are placed in a cubic box with periodic boundary conditions. In our simulations 
we have considered a relatively small number of particles, i.e. iV = 66. Previous studies 
have shown that such a small sample is quite adequate to show interesting off-equilibrium 
behaviour. 

The main quantity on which we will concentrate our attention is the diffusion of the 
particles, i.e. 

A(* , t w ) = ^=i,N\Mt w )-Mt w + t))\ 2 (12) 

The usual diffusion constant is given by D = lim^oo A(t, t w )/t. 

The other quantity which we measure is the response to a force. We add at time t w the 
term ef ■ x&, where / is vector of squared length equal to 3 (in three dimensions) and we 
measure the response 



de 



6=0 



for sufficiently small e. 

The usual fluctuation theorem tells that at equilibrium f3A = R. This relation holds in 
spite of the fact that A(t,t w ) is not the product of two observables one at time t w , the other 
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Figure 1: (3 A (points) and R (crosses) as function the ratio ln(l + s) = ln(l +t/t w ) at V — 1.7 
for t w = 2048. 



at time t w + 1. However it can be written as 

A(t, O = J2 Mt w ) 2 + Mt w + t)f - 2xi(t w ) ■ Xi(t w + t))/N (14) 



i=lJV 



and a detailed analysis JTJJ shows that the fluctuation dissipation theorem is valid also in this 
case. 

In the following we will look for the validity in the low temperature region of the generalized 
relation j3X(A)A = R. This relation can be valid only in the region where the diffusion 
constant D is equal to zero. Strictly speaking also in the glassy region D ^ 0, because 
diffusion may always happens by interchanging two nearby particles (D is different from zero 
also in a crystal); however if the times are not too large the value of D is so small in the glassy 
phase that this process may be neglected in a first approximation. 

We have done simulations for N = 66. For iV = 66 we average over 250 samples in absence 
of the force and on 1000 samples in presence of the force. We have done simulations for 
T = 1.4, 1.5. ..2.0 . The dynamics was implemented using the Monte Carlo method, with an 
acceptance rate fixed around .4. In order to decrease the error on the determination of R we 
follow (as suggested in [0]) the method of computing in the same simulation the response 
function for different particles [0. In other words we add to the Hamiltonian the term 



e J2k=i,N fk ■ x fc where the / are random Gaussian vectors of average squared length equal to 
3. The quantity R can be computed as the average over / of 

iVe • (15) 

The value of e should be sufficient small in order to avoid non linear effects: we have done 
extensive tests for e = .1 and e = .2. We found that e = .2 is in the linear region, but we have 
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Figure 2: R versus /3A at F = 1.6 for t w = t w = 8192. and t w = 2048 

followed the more conservative option e = .1. We present the results for t w = 2048. We have 
done also some simulations at t w = 8192, but we have not observed any systematic shift. 

In fig.l we show the dependence of A and R on the ratio s = t/t w in the low temperature 
region, i.e. at T = 1.7 for t w = 2048. They coincide in the small s region, where FDT holds, 
but they differ for s > .1. 

In fig. 2 we show R versus f3A at t w = 2048, still at T = 1.7. We see two linear regions with 
different slope as follows from the assumption of one step replica symmetry breaking. The 
slope in the first region is 1, as expected form the FDT theorem, while tre slope in the second 
region is definitively smaller that 1. 

Also the data at different temperatures for all values of T > 1.5 show a similar behaviour. 
The value of R in the region where the FDT relation does not hold can be very well fitted by 
a linear function of A as can be seen in fig. 2. The region where a linear fit (with m < 1) is 
quite good correspond to t/t w > .2. The fitted value of m = dR/d((3A) is displayed in fig. 3. 
When m becomes equal to 1, the fluctuation dissipation theorem holds in the whole region and 
this happens at higher temperatures. The straight line is the prediction of the approximation 
m(T) = T/T c , using r c = 1.45 which seems to be rather good.. 

The value of m we find at T — 1.8 (i.e. m = .33 ± .04) is compatible with the value 
(m = .3 ± .1) of ref. extracted from the fluctuation of the stress energy tensor. The 
method described in this note is much more accurate for two reasons: 

• The quantities which we consider becomes larger when we enter in the OFDR region: 
they increase (not decrease) as function of time. 

• The correlation quantity we measure is an intensive quantity which becomes self aver- 
aging in the limit of infinite volume. 

It is amusing to notice (as stressed to me by G. Ruocco) that the simple aging relation 
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Figure 3: The quantity m = as t w = 2048 as function of the temperature. The straight 
line is the prediction of the approximation m(T) = T/T c . 

A(t w ,t w ) = cost for large t w implies that the particle move in average by a constant amount 
in each interval of time 2 K < t < 2 K+1 . If we assume that the movement in each time interval 
are uncorrelated, it follows that in the glassy phase A(t,t w ) oc ln(l + t/t w ) for not too small 
t/t w . This is what happens outside the FDT region, (i.e. t/t w > .2), as it can be seen from 
fig. 1. In fig. 4 we show the data for t w = 1, i.e. the distance from the initial configuration. 
The data seem to display a very nice logarithmic behaviour. 

All these results are in very good agreement with the theoretical expectations based on 
our knowledge extracted from the mean field theory for generalized spin glass models. The 
approximation m(T) = T/T c seems to work with an embarrassing precision. We can conclude 
that the ideas developed for generalized spin glasses have a much wider range of application 
than the models from which they have been extracted. It likely that they reflect quite general 
properties of the phase space and therefore they can be applied in cases which are quite 
different from the original ones. 

The most urgent theoretical task would be now to develop an analytic theory for glasses 
in the low temperature region from which one could compute the function m(t). This goal 
should not be out of reach: a first step in this direction can be found in |H? | 
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Figure 4: The distance from the initial configuration as function of the time in a logarithmic 
scale. 
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